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By detecting redshift drift in the spectra of Lyman-a forest of distant quasars, Sandage-Loeb (SL) 
test directly measures the expansion of the universe, covering the “redshift desert” of 2 < 2 < 5. 
Thus this method is definitely an important supplement to the other geometric measurements and 
will play a crucial role in cosmological constraints. In this paper, we quantify the ability of SL 
test signal by a CODEX-like spectrograph for constraining interacting dark energy. Four typical 
interacting dark energy models are considered: (i) Q — 'yHpc, (ii) Q — 'yHpde, (in) Q ~ 'yHopc, 
and (iv) Q — 'yHopde- The results show that for all the considered interacting dark energy models, 
relative to the current joint SN+BAO+CMB+JTo observations, the constraints on Qrn and Ho 
would be improved by about 60% and 30-40%, while the constraints on w and 7 would be slightly 
improved, with a 30-yr observation of SL test. We also explore the impact of SL test on future 
joint geometric observations. In this analysis, we take the model with Q = 'yHpc as an example, 
and simulate future SN and BAO data based on the space-based project WFIRST. We find that 
in the future geometric constraints, the redshift drift observations would help break the geometric 
degeneracies in a meaningful way, thus the measurement precisions of Qm, Hq, w, and 7 could be 
substantially improved using future probes. 

I. INTRODUCTION 

Since interactions are ubiquitous in nature, it is rather 
natural to imagine that dark energy might directly inter¬ 
act with cold dark matter. Actually, that there is no any 
direct interaction between dark energy and dark matter 
is an additional, strong assumption. Currently, one of 
the most important missions in the field of dark energy 
research is to provide positive/negative evidence for (i.e., 
certify/falsify) the scenario of interacting dark energy in 
the light of observational data. Synthetically using the 
measurements of expansion history and growth of struc¬ 
ture to consistently test the scenario is a fairly important 
way. However, for the scenario of interacting dark energy, 
it is difficult to finely test models using the measurements 
of growth of structure due to some complexness, such as 
the diversity of the construction of covariant 4-vector in¬ 
teraction, the large-scale gravity instability, and the lack 
of abundant, highly accurate data of growth of struc¬ 
ture. Although the progress has been made in this as¬ 
pect since the parametrized post-Friedmann theoretical 
framework for interacting dark energy was proposed and 
applied [1, 2], obstacle due to other factors still exists. 

Under such circumstance, in this work, we only consider 
to use the geometric measurements to constrain the in¬ 
teracting dark energy models; in particular, we focus on 
the future redshift drift data. 

As a purely geometric measurement, Sandage-Loeb 
(SL) test will be crucial to probe the “redshift desert” 

(2 ^ ^ 5) by directly measuring the expansion of 

the universe. It was firstly proposed by Sandage [3] 
to directly measure the variation of redshift of distant 
sources. Then Loeb [4] found a realistic way of detect- 
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ing redshift drift in the spectra of Lyman-a forest of dis¬ 
tant quasars (QSOs). The 39-meter European Extremely 
Large Telescope (E-ELT) in built was equipped with 
a high-resolution spectrograph called CODEX (COsmic 
Dynamics Experiment), which was designed to achieve 
this goal. Cosmological constraints with SL test have 
been studied by numerous works [5-13], some of which 
simulated 240 or 150 quasars to be observed. However, 
as pointed in Ref. [14], only about 30 quasars are bright 
enough or lying at high enough redshift for the actual 
observation, based on a Monte Carlo simulation using a 
telescope with a spectrograph like CODEX. Besides, as 
far as we know, in the most existing papers, the best- 
fit ACDM model to current data is chosen as the hducial 
model, based on which SL test data are simulated. Thus, 
when these SL test data are further combined with other 
actual data to constrain dark energy models, tension be¬ 
tween the simulated SL data and other actual data may 
occur. This is inappropriate and may not give convinc¬ 
ing conclusion on the impact of future SL test data on 
cosmological constraints. 

To avoid inconsistency in data, in our previous 
work [15], we suggested that the best-fit model in study 
to current actual data is chosen as the fiducial model in 
simulating 30 mock SL test data. To give a typical ex¬ 
ample, we only focused on the dark energy model with 
constant w (referred to as rcCDM model). In our recent 
work [16], we extended the discussion to time-evolving 
dark energy model and explored the impact of SL test 
data on dark energy constraints in the future geometric 
measurements. 

Though in Ref. [15], a preliminary SL test analysis has 
been made for two simple interacting dark energy models, 
a synthetical analysis in depth in quantifying the impact 
of future redshift drift data on testing different types of 
interacting dark energy models is still absent. This pa¬ 
per will provide such an analysis. In this paper, we will 



2 


quantify the constraining power of future SL test data on 
different interacting dark energy models, and show how 
the SL test impacts on the parameter estimation. 

For interacting dark energy models, the energy balance 
equations for dark energy and cold dark matter are 

Pde + 'iHpdei^+w) =-Q, (1) 

Pc + 3 iLpc = Q) ( 2 ) 

where pde and pc are the background energy densities 
of dark energy and cold dark matter, respectively. The 
Hubble parameter H = a/a describes the expansion rate 
of the universe and the interacting term Q describes the 
energy transfer rate between dark energy and dark mat¬ 
ter densities. We consider four typical interacting dark 
energy models: (i) Q — jHpc, with which the model is 
called the IwCDMl model, (ii) Q = 'yHpde, with which 
the model is called the IryCDM2 model, (iii) Q = 'yHopc, 
with which the model is called the IicCDMS model, and 
(iv) Q = "fHopde, with which the model is called the 
IwCDM4 model. Here 7 is the dimensionless coupling 
parameter, and the equation-of-state parameter w is con¬ 
sidered to be a constant in this paper. 

In fact, when the CODEX experiment is ready to de¬ 
liver its redshift drift data, other future geometric mea¬ 
surements data will also be available. Therefore, it seems 
that a more meaningful issue is to ask what role the SL 
test will play in parameter estimation for interacting dark 
energy models using the future geometric measurements. 
To address this issue, we also simulate the future SN and 
BAO data based on the long-term space-based project 
WFIRST (Wide-Field Infrared Survey Telescope). It is 
quite common that program in built is changed or post¬ 
poned. Besides, because science and technology develop 
rapidly, it is hard to well and truly quantify the per¬ 
centages of the parameter estimation by SL test data. 
We only take WFIRST as an example. Owing to the 
fact that it covers the redshift desert (2 < z < 5), SL 
test is a valuable supplementary to other geometric mea¬ 
surements. Moreover, because of the different degener¬ 
acy orientations of SL test and other geometric measure¬ 
ments, it is very credible that SL test can effectively break 
the strong degeneracy in other geometric measurements 
and greatly improve the corresponding cosmological con¬ 
straint results. It’s necessary to point out that this anal¬ 
ysis does not mean the necessity of combining the data 
from WFIRST and CODEX in the future, and that we 
merely explore the ability of redshift drift to break degen¬ 
eracy in other geometric measurements and to improve 
the accuracy of cosmological constraints. 


II. METHODOLOGY 

Our procedure is as follows. Interacting dark energy 
models are first constrained by using the current joint 
SN-|-BAO-|-CMB-fiLo data, and then for each case the 


best-fit model is chosen to be the hducial model in pro¬ 
ducing the simulated mock SL test data. The obtained 
SL test data are thus well consistent with the current 
data. Therefore, it is rather appropriate to combine the 
mock SL test data with the current data for further con¬ 
straining the interacting dark energy models. We per¬ 
form an MCMC likelihood analysis [17] to obtain 0(10®) 
samples for each model. 

For current data, the most typical geometric measure¬ 
ments are chosen, i.e., the observations of SN, BAO, 
CMB, and Hq. For the SN data, the SNLS compila¬ 
tion [18] with a sample of 472 SNe is used. For the 
BAO data, we consider the rsfDv{z) measurements from 
6 dFGS {z = 0.1) [19], SDSS-DR7 {z = 0.35) [20], SDSS- 
DR9 {z = 0.57) [21], and WiggleZ {z = 0.44, 0.60, and 
0.73) [22] surveys. For the CMB data, we use the Planck 
distance posterior given by Ref. [23] . As dark energy only 
affects the CMB through the comoving angular diameter 
distance to the decoupling epoch (and the late-time ISW 
effect), the distance information given by the CMB dis¬ 
tance posterior is sufficient for the joint geometric con¬ 
straint on dark energy. We also use the direct measure¬ 
ment result of the Hubble constant in the light of the 
cosmic distance ladder from the HST, Hq = 73.8 ± 2.4 
km s“^ Mpc“^ [24]. 

Next, we briefly review how to simulate the SL test 
data. This method is just to directly measure the red¬ 
shift variation of quasar Lyman-a absorption lines. The 
redshift variation is defined as a spectroscopic velocity 
shift [4], 


Av = 


Az 
1 + z 


— HgAto 



E{z) 
1 + z 


( 3 ) 


where Ato is the time interval of observation, and E{z) = 
H{z)/Hq is given by specihc dark energy models. 

According to the Monte Carlo simulations, the un¬ 
certainty of Av expected by CODEX can be expressed 
as [14] 


= 1-35 


S/NY^ 

2370) 


Nqso 

30 


- 1/2 


1 + ^QSO 

5 


X 

cm s 


( 4 ) 

where S/N is the signal-to-noise ratio defined per 0.0125 
A pixel, Nqso is the number of observed quasars, zqso 
represents their redshift, and the last exponent x = —1.7 
for 2 < z < 4 and x = —0.9 for z > 4. 

To simulate the SL test data, we first constrain the in¬ 
teracting dark energy models by using the current data 
combination. The obtained best-ht parameters are sub¬ 
stituted into Eq. (3) to get the central values of the SL 
test data. We choose Nqsq = 30 mock SL data uniformly 
distributed among six redshift bins of z:qso G [2,5] and 
typically take Ato = 30 yr in our analysis. The error 
bars are computed from Eq. (4) with S/N = 3000. 

In order to quantify the power of SL test in future 
high-precision joint geometric constraints on dark energy, 
we simulate future SN and BAO data based on the long¬ 
term space-based project WFIRST using the method pre- 
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sented in Ref. [25], and take the interacting dark energy 
model with Q = jllpc as an example. We simulate 2000 
future SNe distributed in 16 bins over the range z = 0.1 
to z = 1.7. The observables are apparent magnitudes 
rrii = M + p{zi), where M represents the absolute magni¬ 
tude, and p,{zi) is the distance modulus. We also include 
an additional “near sample” of 500 SNe at z « 0.025. For 
future BAO data, we simulate 10000 mock BAO data uni¬ 
formly distributed among 10 redshift bins of z S [0.5, 2], 
with each Az^ centered on the grid z^. The observables 
are expansion rate H{z) and comoving angular diameter 
distance d‘^{z) = dL{z)/{l + z). For details, we also refer 
the reader to Ref. [16]. 


III. RESULTS AND DISCUSSION 


First, we constrain the icCDM model and four typ¬ 
ical interacting dark energy models from the current 
SN-|-BAO-|-CMB-|-i7o data combination, and present the 
detailed fit results in Table I. From this table, one can 
clearly see that for wCDM model, u> < — 1 is preferred 
at about the l.Scr level, while w < — 1 is preferred at 
more than 2.2(t level for all the four interacting dark en¬ 
ergy models. For the wCDM, IuiCDM2, IwCDMS and 
IwCDM4 models, ilch^ can be tightly constrained, and a 
smaller value is preferred. But for the IwCDMl model, 
cannot be well constrained, and a bigger value is 
more favored in this case. The coupling 7 is tightly 
constrained in the IwCDMl model, but its constraint 
is much weaker in the IwCDM2, IwCDMS and IwCDMd 
models. For IwCDMl model, 7 < 0 is preferred at about 
2.1cr level, while 7 < 0 is slightly favored at about 1 .do- 
level for IwCDM2, IwCDMS and IwCDM4 models. 

Then we combine the simulated 30-yr SL test data 
with current data and show the constraint results in 
Figs. 1-3. To give an intuitive comparison, the re¬ 
sults from current only data are also presented. Fig¬ 
ures 1 and 2 show the joint constraints on the u>CDM, 
IwCDMl, It(;CDM2, IwCDM3 and IwCDM4 models in 
the rtm-Ho and flm-w planes, respectively. Figure 3 
shows the joint constraints on the four interacting dark 
energy models in the plane. In these three hgures, 

the 68.3% and 95.4% CL posterior distribution contours 
are shown. The data combinations used are the current 
only and the current-l-SL 30-yr combinations, and their 
constraint results are shown with red and blue contours, 
respectively. Here we use “current” to denote the current 
SN-|-BAO-|-CMB-fi7o data combination for convenience. 
One can clearly see that the degeneracies are well broken 
with the SL test data for all the models. 

The 1(7 errors of the parameters w, 7 , 0^, and Hq for 
the five models with current only and current-l-SL 30-yr 
data are given in Table 11. With the 30-yr SL observa¬ 
tion, the constraints on 0 ^ and Hq will be improved, 
respectively, by 81.4% and 64.6% for the rcCDM model, 
by 68.3% and 27.9% for the IwCDMl model, by 58.0% 
and 44.1% for the Ir(;CDM2 model, by 58.3% and 44.7% 







FIG. 1: Constraints (68.3% and 95.4% CL) in the Qm- 
Ho plane for wCDM, IwCDMl, I«)CDM2, IwCDM3, and 
IwCDM4 models with current only and current-|-SL 30-yr 
data. 


for the Iu>CDM3 model, and by 58.3% and 43.3% for the 
ItcCDM4 model. Therefore, we can see that with a 30- 
yr observation of the SL test the geometric constraints 
on dark energy would be improved enormously. For all 
the considered interacting dark energy models, the con¬ 
straints on Qm and Hq would be improved, relative to 
the current joint observations, by about 60% and 30- 
40%, with the SL 30-yr observation. 

We also discuss the impact of the SL test data on con¬ 
straining the equation of state w and the coupling 7 . 
The SL 30-yr observation helps improve the constraints 
on w by 24.4%, 20.5%, 14.6%, 10.9% and 13.1% for 
the wCDM, ItcCDMl, IwCDM2, IwCDM3 and IwCDM4 
models, respectively. The SL 30-yr observation helps im¬ 
prove the constraints on 7 by 30.5%, 9.0%, 9.7% and 
8 . 8 % for the IwCDMl, Ii(;CDM2, IwCDM3 and IwCDM4 
models, respectively. Therefore, we find that among the 
four interacting dark energy models, the IwCDMl model 
is the best one in the sense that the constraint results 
could be improved by the SL test data. 

We have discussed the quantification of the impact 
of future SL test data on constraining interacting dark 
energy from current SN-l-BAO-|-CMB-fiJo observations. 
The results show that future SL test data can effectively 
break the degeneracies in the current data for interacting 
dark energy models and thus will provide fairly important 
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TABLE I: Fit results for the uiCDM, IwCDMl, IwCDM2, IwCDMS and IwCDM4 models using the current 
CMB+BAO+SN+ii'o data. 


Parameter 

wCDM 

IwCDMl 

Im;CDM2 

IwCDMO 

I?nCDM4 

Q.hh? 

Q.ch? 

w 

7 

ffo 

n nooi 

U.UZZ 10 _Q 00029 

U.1ZU1_Q QQ22 

-1.103 + 0.058 

0.2844t°;°J°^3 

70.7411:26 

A AAAAQ + O.00033 

u.uzzzy_o.00024 

A 1 9OQ+0.0070 

U.iZ50_o.o049 

1 1 0/5+0.062 
l.lCiD_o.055 

-0.01121^0?! 

0.2849l°:°“l 

72.8111:61 

0.02232 ± 0.00027 

u.r^U3_Q QQ22 

ir9+0.064 
l.iOZ_o 072 

-0.0284lHll,®6 

0.2834lg:i;?l4 

71.0911:1^ 

0.02231 ± 0.00028 

A 1 91 r+0.0020 
U.lZlO_o QQ23 

-i-isoIHm 

-0.0280lg:°«l 

A 9O9r+0.0096 

U.ZOZO_o 0100 

71.34 ± 1.33 

A A 990 A+ 0 .00030 

U.UZZOU_o 00026 

U.1ZUO_0 0019 

i.i4y_0.077 

—0 

u.uzi:;t;_0.0264 

n 9099 + 0.0102 
U.ZOZZ_o 0095 

71.1511:^^ 


TABLE II: Errors of parameters in the wCDM, IiuCDMl, IiuCDM2, ItnCDMS and IwCDMd models for the fits to current only 
and current+SL 30-yr data. 


current only 

current + SL 30-yr 

Error wCDM IwCDMl IwCDM2 IwCDMO IwCDM4 

wCDM IwCDMl Im;CDM2 IwCDM3 IwCDM4 


a(w) 

0.082 

0.083 

0.096 

0.092 

0.099 

0.062 

0.066 

0.082 

0.082 

0.086 

< 7 ( 7 ) 

- 

0.0095 

0.0300 

0.0279 

0.0340 

- 

0.0066 

0.0273 

0.0252 

0.0310 


0.0140 

0.0142 

0.0143 

0.0139 

0.0139 

0.0026 

0.0045 

0.0060 

0.0058 

0.0058 

<r(i7o) 

1.81 

2.51 

1.77 

1.88 

1.78 

0.64 

1.81 

0.99 

1.04 

1.01 




current only 


current + SL 30-yr 



FIG. 3: Constraints (68.3% and 95.4% CL) in the Q.m-'y plane 
for ItnCDMl, IwCDM2, ItnCDM3, and IwCDMd models with 
current only and current+SL 30-yr data. 


FIG. 2: Constraints (68.3% and 95.4% CL) in the Q,m-w plane supplement to the other observations. In the following, 
for wCDM, IiuCDMl, IwCDM2, IwCDM3, and IwCDM4 we will further explore what role the SL test will play 
models with current only and current+SL 30-yr data. in the future geometric constraints on interacting dark 

energy. We take the IwCDMl model (with Q = ^Hpc) 
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wCDM future only future + SL 30-yr 



FIG. 4: Constraints (68.3% and 95.4% CL) in the Qm-Ho 
plane and in the Qm-w plane for tnCDM model with future 
only and future+SL 30-yr data. 


IwCDMI (Q= 7 HpJ future only future + SL 30-yr 


72.90 

72.86 


-1.05 

-1.09 

V 

1° 

72.82 


-1.13 

X 

72.78 


-1.17 

X 

72.74 


-1.21 



0.26 0.28 0.30 0.32 ' 0.26 0.28 0.30 0.32 




FIG. 5: Constraints (68.3% and 95.4% CL) results for 
IwCDMI model with future only and future-fSL 30-yr data. 


as an example for this analysis and compare the result 
with that of the wCDM model. As mentioned above, we 
simulate future SN and BAO data based on the WFIRST 
mission in this analysis. 

Figure 4 shows the joint constraints on the wCDM 
model in the flm-Ho and flm-w planes. Figure 5 shows 
the joint constraints on the IwCDMI model in the ilm- 
Hq, fljn-w, and w-y planes, respectively. The 

68.3% and 95.4% CL posterior distribution contours are 
shown. The data combinations used are the future only 
and the future -I- SL 30-yr combinations, and their con¬ 
straint results are shown with red and blue contours, re¬ 
spectively. The la errors of the parameters w, 7 , fim, 
and ffo for the wCDM model and the IwCDMI model for 
the above data combinations are given in Table III. Note 
that here we use “future” to denote the data combina¬ 
tion of future SN and BAO for convenience. It is shown 


TABLE III: Errors of parameters in the wCDM and IwCDMI 
models for the fits to future only and future-|-SL 30-yr data. 


Error 

future only 

future -|- SL 30-yr 

wCDM IwCDMI 

wCDM 

IwCDMI 

a{w) 

0.0083 

0.0416 

0.0067 

0.0132 

cr( 7 ) 

- 

0.1119 

- 

0.0181 


0.0021 

0.0190 

0.0016 

0.0045 

a{Ho) 

0.0296 

0.0374 

0.0286 

0.0315 


that with the 30-yr SL observation, the constraints on 
fim and ffo will be improved by 23.8% and 3.4% for the 
wCDM model, and by 76.3% and 15.8% for the IwCDMI 
model. For the wCDM model, the constraints on w can 
be improved by 19.3%, with the SL 30-yr observation, 
while for the IwCDMI model, the SL 30-yr observation 
helps improve the constraints on w and 7 by 68.3% and 
83.8%, respectively. 

Comparing Figs. 4 and 5, we find that in the fu¬ 
ture geometric constraints, the redshift drift observa¬ 
tion could not break the degeneracies for the wCDM 
model, but could efficiently break the degeneracies for the 
IwCDMI model. In the future geometric constraints, for 
the IwCDMI model, the SL 30-yr observation would help 
improve the measurement precisions of 11 ^, -ffoi ^ and 7 
by more than 75%, 15%, 65%, and 80%, respectively. 

The main purpose of this work is to investigate the 
possible impact of future SL test data on existing ge¬ 
ometric measurements. In the future, in case ACDM 
model is excluded, one important work is to constrain 
the coupling parameter 7 , if we wish to find evidence for 
the existence of interaction between dark sectors. There¬ 
fore, it is quite meaningful to investigate the impact of 
future SL test data on parameter estimation for interact¬ 
ing dark energy models. For all considered interacting 
wCDM models in this work, it is shown that SL test can 
effectively break the existing parameter degeneracies and 
greatly improve the precisions of parameter estimation. 
The results are consistent with the cases of the ACDM, 
wCDM, and wqWqCDM models [15, 16]. By considering 
more models, we can conclude that the improvement of 
parameter estimation by SL test data is independent of 
the cosmological models in the background. This shows 
the importance of including SL test data in future cos¬ 
mological constraints. 


IV. SUMMARY 

In this paper, we have analyzed how the redshift drift 
measurement (i.e. SL test signal) would impact on pa¬ 
rameter estimation for the interacting dark energy mod¬ 
els. By detecting redshift drift in the spectra of Lyman-a 
forest of distant quasars, SL test directly measures the 
expansion rate of the universe in the redshift range of 
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2 < z < 5, providing an important supplement to other 
probes in dark energy constraints. We consider four typ¬ 
ical interacting dark energy models: (i) Q = yiJpcj (ii) 
Q = -fHpde, (iii) Q = 'yHoPc, and (iv) Q = jHopde- 

Following our previous works [15, 16], in order to guar¬ 
antee that the simulated SL test data are consistent with 
the other geometric measurement data, we used the best¬ 
fitting dark energy models constrained by the current 
combined geometric measurement data as the hducial 
models to produce the mock SL test data. 

We showed that the SL test data are extremely helpful 
in breaking the existing parameter degeneracies. Com¬ 
pared to the current SN-|-BAO-|-CMB-fiJo constraint re¬ 
sults, the 30-yr observation of SL test could improve the 
constraints for all the considered interacting dark energy 
models on and Hq by about 60% and 30-40%, while 
the constraints on w and 7 can be only slightly improved. 

We also quantified the impact of SL test data on in¬ 
teracting dark energy constraints in the future geometric 
measurements. To do this analysis, we simulated the 


future SN and BAO data based on the long-term space- 
based project WFIRST. We found that the SL test could 
also play a crucial role in the future joint geometric con¬ 
straints, especially for the constraints on w and 7 . Taking 
the interacting dark energy model with Q — "fHpc as an 
example, the 30-yr observation of SL test would help im¬ 
prove the measurement precision of Hq, w and 7 by 
more than 75%, 15%, 65%, and 80%, respectively. 
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